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We consider fermion-dimer scattering in the presence of a large positive scattering length in the 
frame of functional renormalization group equations. A flow equation for the momentum dependent 
fermion-dimer scattering amplitude is derived from first principles in a systematic vertex expansion of 
the exact flow equation for the effective action. The resummation obtained from the nonperturbative 
flow is shown to be equivalent to the one performed by the integral equation by Skorniakov and 
Ter-Martirosian (STM). The flow equation approach allows to integrate out fermions and bosons 
simultaneously, in line with the fact that the bosons are not fundamental but build up gradually 
as fluctuation induced bound states of fermions. In particular, the STM result for atom-dimer 
scattering is obtained by choosing the relative cutoff scales of fermions and bosons such that the 
fermion fluctuations are integrated out already at the initial stage of the RG evolution. 
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The scattering of three nonidentical fermions has been considered long time ago by Skorniakov and Ter-Matrirosian 
(STM) l]. For short-ranged interactions and a positive scattering length, the two-body sector supports a shallow 
bosonic bound state, the dimer. If the scattering length is large compared to the effective range, then the low energy 
physics becomes universal in the sense that all observables may be expressed in terms of a single length scale, the 
fermionic scattering length a. Making use of the insensitivity with respect to the short distance physics, STM extracted 
the universal ratio of fermion-dimer to fermion scattering length to be a^/a = 1.18. 

The three-body problem has a long history. In the 70s, Efimov extended the results for fermions to the bosonic 
case, where a short distance scale is needed to stabilize the particles, giving in turn rise to the Efimov effect, i.e. 
the existence of a sequence of three-body bound states Q. In subsequent papers, he revealed the universality of the 
three-body problem in systems with large scattering length [3| , and furthermore considered effective range corrections 

- n 

The first derivation of the STM equation using Feynman diagrams was performed in [5[ using a purely fermionic 
formulation. In the 90s, the three-body problem was considered in the light of Effective Field Theory. The problem 
was reformulated in terms of an effective theory, where the dimer degree of freedom is implemented explicitly 0, 0] ; 
for a review see Extensions of the standard STM problem in this framework include higher order corrections 
in the effective range Q and higher partial waves Recently, the problem has seen renewed interest, developing 
into an important nontrivial benchmark for new techniques. Those include nonperturbative methods like quantum 
mechanical approaches in position space [101 ] , as well as perturbative techniques like the e expansion around the critical 
dimension (ill. 




In this paper, we address the problem of fermion-dimer scattering in the frame of nonperturbative functional 
renormalization group equations (FRG) for the effective action 0jE1I3- This technique has been used successfully 
to quantitatively investigate critical phenomena, i.e. to analyze the universal physics at very long distances close to a 
phase transition . It has also been applied for the study of complex many-body systems [r| [It], , where typical 
length scales are set by the mean interparticle spacing and the inverse temperature. In this paper we demonstrate how 
it can be used to address nonperturbative fluctuation problems on even much shorter distances, i.e. the scattering 
of few particles. Our approach is based on a systematic vertex expansion of the effective action, which keeps the full 
momentum dependence of the one-particle irreducible (1PI) vertices. 

Vertex expansion schemes have been worked out in various applications [l||. In these works, the resulting flow 
equations exhibit a complexity which necessitates a full numerical treatment. Here we take advantage of the fact 
that for the few-body scattering problem at low energies, the complicated momentum and frequency dependence 
can be considerably simplified. This also allows for a direct comparison of our method to other schemes. Focusing 
on three-body (fermion-dimer) scattering here, we stress that the method is also suitable to address four-particle 
(dimer-dimer) scattering (20l. [2lj. which will be discussed in a future publication. Furthermore, we point out that the 
diagrams resummed during the flow involve both inner fermion and boson lines. Situations like this also appear in 
various man y-b ody problems, e.g. in the analysis of quantum critical points where both fermions and bosons develop 
zero modes [22|]. The STM problem allows to study the simultaneous elimination of fermion and boson degrees of 
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freedom in a relatively controlled setting. It may thus provide valuable hints for the future treatment of coupled 
fermion-boson theories in the frame of functional RG, as initiated in recent work [23j | . 

This paper is organized as follows: In Sect. m we briefly sketch the FRG method and formulate the problem in 
terms of a two-channel model for both a stable fermion field and a composite boson field, coupled via the Feshbach 
coupling h v . In the limit of large Feshbach couplings, universality emerges naturally as an infrared stable fixed point 
for the renormalized Feshbach coupling [13, HH > and our model becomes equivalent to a single channel model with 
point-like fermionic interactions. At this stage the problem is formulated for arbitrary density and temperature. In 
Sect. IIII1 we then specify the prescription to project onto the physical vacuum of vanishing density and temperature, 
as appropriate for few-body scattering. The following section is devoted to a discussion of the (exact) solution of 
the two-body problem in the FRG framework. Sect. IVlthen deals with fermion-dimer scattering: the FRG equation 
for the fermion-dimer vertex is derived and cast into a form closely resembling the STM integral equation. We then 
consider solutions of the RG equations for various choices of the relative cutoff scale at which fermions and dimers 
are integrated out relatively to each other, and show in which limit the STM result is reproduced. Conclusions are 
drawn in Sect. |VU 



II. METHOD AND APPROXIMATION SCHEME 

We study the scale dependence of the effective average action Tk for reviews see [l3|, EH- It includes all 
fluctuations with momenta q 2 > k 2 . In the limit k — ► where the averaging scale k is removed, all fluctuations 
are included and approaches the full effective action. In practice, the scale dependence is implemented by 

introducing suitable cutoff functions Rk{q) in the inverse propagators. The dependence of T k on k obeys an exact 
flow equation, 

d k T k = isTr(ri 2) +R k )- 1 d k R k = ~STrd k log(rW + R k ). (1) 

Here, the "supertrace" STr sums over spatial momenta q and Matsubara frequencies to as well as over internal indices 
and species of fields, with a minus sign for fermions. The effective action is formulated in Euclidean spacetime. The 

(2) . . (2) 

second functional derivative T k represents the full inverse propagator in the presence of the scale k. Both T k and T k 
are functional of the fields. For the last equation, the derivative d k is defined to act on the explicit scale dependence 
set by the cutoff function R k , and not on the implicit scale dependence of T^. This leads to a compact notation, and is 
advantageous to make direct contact with diagrammatic representations. In practice the above functional differential 
equation can only be solved approximately using a suitable truncation of the full effective action functional. The 
vertex expansion is obtained from expanding the last expression in powers of the fields. 

The effective action is the generating functional of the 1PI n-point correlation functions. In the physical vacuum 
state of vanishing density and temperature, these objects can be directly related to the scattering amplitudes. 

The effective action for fermions interacting via a Feshbach resonance can be described by a simple two-channel 
ansatz. In momentum space, and after analytical continuation to Euclidean frequencies (wm — > —iuj, where wm is the 
Minkowski frequency) it reads 

a A )i)(Q) + <p*(Q)P v (QMQ)] (2) 

Q 2 - Q3)(<P*(Qi)MQ2)MQa) - p(QiWi(W2(Q3)) 

with four-momentum Q = (uj,q). Here if) = (^1,^2) represents the stable nonrelativistic fermionic atom field. We 
further introduce a composite boson field tp which mediates the interactions between the fermions via the Yukawa or 
Feshbach coupling h m. T he composite field can play various physical roles, depending on the region of parameter space 
under consideration [T7| . as well as on the averaging scale k. Here we will restrict ourselves to the vacuum limit where 
the scales associated to many-body physics vanish (T = n = 0), as well as to the limit of broad Feshbach resonances, 
where the value of the (bare) Feshbach coupling drops out as an independent scale in the problem, h v — ► 00. In this 
limit, the bosonic field is purely auxiliary at a high scale k with no propagation, and acquires dynamics only in the 
limit k — > by the virtue of fluctuations, as will be explained in more detail below. The only remaining scale in 
the problem is then the fermionic scattering length a. This realizes "large" scattering lengths, in the sense that the 
observable physics is uniquely determined by the latter value, becoming insensitive with respect to further microscopic 
information. This has an important aspect of universality, which is discussed for the few body-problem in [8| , and 
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with an emphasis on the implications for the many-body problem in [2J, [25|, [26[ . The projection procedures on the 
physical vacuum and on the universal broad resonance limit are specified in the next section. 

Our units are h = k^ = c = 1. Furthermore, we measure all momenta in units of some reference scale k, and energies 
in units of i k — k /2M, where M is the nonrelativistic mass of the fermions. This leads to a dimensionless scaling 
formulation of the effective action [§f|. We have the following relations between dimensionless and dimensionful 
(denoted with a hat) quantities , 

Q = {uj,q) = (w/e k ,q/k), a A = a A /e k , h lp = 2Mk-^%, P 9 {Q) = P v (fi) / 'e k , V = k' 3 / 2 ^, cp = k^' 2 <p. (3) 

The physical meaning of the parameter a A again depends on the parameter regime under consideration. In the 
many-body context at finite density and temperature, it plays the role of the chemical potential for the fermions. In 
the vacuum, it represents half the binding energy of a dimer (cf. Sect. HVj) . which is nonzero for positive scattering 
lengths. 

The inverse bosonic propagator consists of two contributions: a classical or "bare" part v and a fluctuation part 
SP V , 

P V (Q) = v + 5P v (Q). (4) 

The parameter v defines the initial condition for the flow of the boson propagator. It includes the physical detuning 
v{B) from the Feshbach resonance as well as a counter term 8v ln needed for the ultraviolet (UV) renormalization of 
the problem as discussed in Sect. [iVj 

v = v{B)+8v ia , v = Z/e k , v{B)^ii{B-B ). (5) 

v{B) measures the distance from the Feshbach resonance at magnetic field Bo, with /i being the effective magnetic 
moment of the atoms in the open channel. Dimensionless detuning and dimensionless scattering lengths (in the 
absence of an open channel background scattering length) are related by 



h 2 



ak. (6) 



8tti/(.B) ' 
The inverse fermion propagator 

Pf(Q) = iLU + q 2 -a A (7) 

does not receive any renormalization corrections in vacuum as will be shown in Sect. IIIII Therefore, we do not 
introduce scale dependent running couplings in the fermionic part of Eq. ^ . The running couplings for this part of 
the effective action are thus h v and P V (Q). 

In order to describe scattering processes involving more than two fermions, we need to extend the truncation Eq. 
@. In particular, the scattering amplitude of a fermion off a dimer is described by the amputated connected part of 
the Green function (0\ipip(p*^\0) 0. Thus we need to include a fermion-dimer coupling 

S(Qi + Q2-Q3- Qi)SX 3 (Qu Q2, Q 3 MQiMQ2)^(Q 3 )4>HQi)- (8) 

The fermion fields are contracted as ipip^ = ipa^aptp^ in spin space, with S a p the identity matrix in two dimensions. 
The order of ip and ip> is important due to the Grassmann nature of the fermionic fields and chosen such that it matches 
the standard conventions for the fermion-dimer scattering amplitude. The coupling depends on three independent 
four-momenta by momentum conservation. Still the four-momentum dependence is very involved and will be largely 
simplified below. However, a point-like truncation of the interaction vertices (no momentum dependence) turns out 
to be insufficient to reach satisfactory precision in this problem. The dimensionless and dimensionful fermion-dimer 
couplings are related by 

SX 3 = 2MkSX 3 . (9) 



In an earlier publication 12611 , the dimensionless quantities were denoted with a tilde superscript, while the quantities without superscript 
were reserved for "renormalized" (rescaled with the wave function renormalization) quantities. We do not introduce such a wave function 
renormalization here, and it is advantageous to reserve the simplest notation for the dimensionless quantities. 
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At this point we stress the systematic nature of the truncation advocated here. The vertex expansion is an expansion 
in the number of the fields. Our truncation is complete up to third order in the fields. At fourth order in the 
fields, there are two more terms which are compatible with U(l) symmetry, namely 

Q 4 )<^MQi, Q 2 , QdtfiQiWiQiWHQaMQ*), (10) 

Qi)5\ A (Qi, Q 2 , Q 3 )<P*(QiMQ2)<P*(Q3MQa). 

The first one describes the scattering of two fermions, the second one the scattering of two dimers, i.e. four-fermion 
scattering in the presence of a bound state. The flow of these couplings has to be taken into account in a systematic 
expansion to fourth order in the fields. However, we will show in Sect. Mil (i) The first term is not generated by the 
flow in the vacuum limit considered here. This means that the fermionic two-body sector is fully described by our 
Feshbach model, (ii) The second vertex describes interactions in the four-body sector and is generated in the vacuum 
limit. In principle, it could couple into the flow of the other vertices. We find that this is not the case in the vacuum 
limit. 

Finally, we specify the regulator functions Rk- We work with a momentum independent, mass- like cutoff function 
for fermions and bosons 

R k , F = k 2 , R kt<p = ck 2 , (11) 

similar to 25]. The choice of the dimensionless number c, which we specify below, sets the relative scale for the 
elimination of the fermionic and bosonic degrees of freedom in the renormalization group flow. The optimal choice of 
c ensures an equal effective cutoff scale for fermions and bosons, where the effective cutoff scale is composed of the 
cutoff function piece plus possible physical mass terms [U [H, 0] ■ 

A scheme mass-like cutoff function is possible if the fermionic "chemical potential" <ta < 0, which is indeed the case 
in the vacuum on the BEC side. For our purpose the mass-like cutoff is advantageous since it allows for most direct 
comparison with conventional diagrammatic techniques. However, for high accuracy calculation an optimized cutoff 
[27l | would be more appropriate. 




5(Qi -Q2 + Q3- 
S(Qi -Q2 + Q3- 



Qi,...Q* 



III. VACUUM LIMIT 



In this paper we consider a specific regime in parameter space where the effective action T — T^o describes the 
scattering of particles in vacuum, which interact via a positive s-wave scattering length. The vacuum projection of 
the effective action is obtained from V in the limit n — > 0, T — ► 0. 

The prescription, which projects the effective action on the vacuum limit reads [26| : 



^vac — lim Tfc^o 



(12) 

T>T C 

Here kp = (37r 2 n) 1//3 is directly related to the density of the system by definition, such that it can be viewed as the 
inverse mean interparticle spacing fcp ~ 1 / d. Taking the limit k-p — > then corresponds to a diluting procedure where 
the density of the system becomes arbitrarily low. However, the limit is constrained by keeping the dimensionless 
temperature T = T/e-p — 2MT/k F above criticality. This ensures that many-body effects such as condensation 
phenomena are absent. Of course, the system becomes arbitrarily cold, since the absolute temperature scales as 
f oc fcf -> 0. 

We find that for n = T = the crossover at finite density turns into a second-order phase transition in vacuum 
[24I I2H ] as a function of the magnetic field B. In order to see this, we consider the momentum independent parts 
in both the fermion and the boson propagator, — a a (the "chemical potential" for the fermions in vacuum) and m 2 ^ , 
which act as gaps for the propagation of fermions and bosons. Here the bosonic mass term is defined as the zero 
frequency and momentum component of the boson propagator, 

ml = P fc =o, v (Q = 0). (13) 

Taking the vacuum limit in the above mentioned form, we find the following constraints, which separate the two 
qualitatively different branches of the physical vacuum [261 ] , 

m 2 > 0, (Ja = atom phase (a -1 < 0), 

m 2 =0, (Ja < molecule phase (a -1 > 0), (14) 
m =0, a a = resonance (a -1 = 0). 
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These formulae have a simple interpretation: On the BCS side, the bosons experience a gap > and the low-density 
limit describes only fermionic atoms. On the BEC side, the situation is reversed: fermion propagation is suppressed 
by a gap — a a- The ground state is a stable molecule, and the fermionic chemical potential can be interpreted as half 
the binding energy of a molecule, cm = 2<ta (26j . — cta is the amount of energy that must be given to each fermion in 
a dimer to reach the fermionic scattering threshold. 

Evaluating the flow equation for the bosonic mass term m 2 , with the constraint (fT4|) on the BEC side (see below), 
one finds the well-known universal relation between binding energy and scattering length in vacuum, cm = — l/(Ma 2 ) 
in the broad resonance limit h v — ► oo. This establishes the second order nature of the vacuum phase transition - 
the resonance at a -1 = is smoothly approached. For finite h Vl scaling violations 0{eyi/hjt > ) emerge 2 . This gives a 
glance at the status of universality related to the value of h v . 

On the technical side, the procedure specified above leads to a massive simplification of the diagrammatic structure 
as compared to the finite density and temperature system. By the aid of the residue theorem, it is straightforward 
to prove the following statement (25^ : All diagrams whose inner lines point in the same direction (thereby forming a 
closed tour) do not contribute to the flow in vacuum. Such diagrams have all poles in the same half of the complex 
plane. The argument holds for frequency and momentum dependent vertices provided that possible poles lie in the 
same half of the complex plane as those of the propagators. We can now analyze the one-loop diagrams which would 
possibly generate an RG flow of the couplings under consideration. An analysis of the one-loop diagrams is sufficient 
due to the one-loop structure of the exact flow Eq. (JTJ. Applying the above statement, we find: 

(i) The fermion propagator is not renormalized in vacuum as can be seen on diagrammatic grounds, and using the 
above argument. 

(ii) A four-fermion vertex ~ (ip^ip) 2 is not generated by the flow. In our model, we have eliminated such a vertex 
to describe fermion-fermion scattering in favor of the coupling to the auxiliary boson degree of freedom (Hubbard- 
Stratonovich transformation). The fact that the vertex is not regenerated by the flow indicates that the Hubbard- 
Stratonovich transformation is very efficient here. We note that at finite density and temperature, such a vertex is 
indeed generated, describing the effect of particle-hole fluctuations. 

(iii) The four-boson (dimer-dimer) vertex Eq. (jTUJ) does not couple into the flow of the couplings we are considering 
in the frame of our truncation in vacuum. This implies that the four-particle sector of the theory does not affect 
the three-particle sector, which is physically sound. On the other hand, we find that the flow in the four-particle 
sector, described by the dimer-dimer vertex, is affected by the fermion-dimer vertex. The same pattern is observed 
for the mutual influence of the two- and three particle sectors. Our vertex expansion therefore seems to respect the 
hierarchy which is expected from physical intuition. At finite density, the interpretation of the vertices as representing 
scattering in sectors with definite particle number is spoiled. As expected, the many-body analogs of these vertices 
then do not respect the hierarchy any more. 



IV. TWO-BODY SECTOR: UV RENORMALIZATION AND UNIVERSALITY 



In order to make contact with experiment, we have to relate the microscopic or bare parameters which characterize 
the theory at a high momentum scale fcj n to the observables for two-atom scattering in vacuum, like the scattering 
length a, the molecular binding energy or an effective range. We therefore choose the initial parameters at the UV 
scale k ln such that the two-body observables are matched in the limit k — > 0. 

In the two-body sector defined by Eq. ^ we need to consider the flow of the Feshbach coupling h v and the inverse 
boson propagator Pk, v (Q)- 

We consider the flow of the Yukawa coupling first. We find 

d k h v = 0. (15) 

The non-renormalization of the Feshbach coupling in vacuum can be traced back to the U(l) symmetry, and extends to 
the case of a momentum dependent Feshbach coupling. This statement holds in the absence of a fermionic background 
coupling. Extending the truncation to take such a coupling into account leads to a renormalization of h v [25l r2t| , 
which is compatible with charge conservation. 

The boson propagator is more involved and we discuss it in detail. The prescription projecting onto this object is 



2 The situation is further complicated in the presence of an additional scale set by a finite background scattering length |2(I . 
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given by the functional derivative 



V k , v {Q 1 ,Q 2 ) 



5 V *{Q 1 )8 V {Q 2 ) 



(16) 



v= 0,ip=0 



We extract it by applying this prescription to both sides of Eq. (Q]) , where in practice we expand the logarithm on 
the rhs in powers of the bosonic field. The inverse bosonic propagator is diagonal in momentum space, Vk^iQ, K) = 
Pk, v {Q)5{Q - K). This yields the flow equation (cf. Fig. QJ 



kd k 6P k , v (Q) = - kd-, 



h 2 



K 



(Pf(-K) + fc 2 ) (Pf(K + Q) + fc 2 ) 



8tt y/iuj/2 + q 2 /4:-<T A + k 2 



(17) 



The flow is large (~ k) for large cutoffs fc. This reflects the presence of a relevant parameter and indicates the 
necessity of an UV renormalization in the language of the flow equation. To make the physics more transparent, we 
may consider the flows of the Q ^ and Q = (mass term) components separately, 



kd k [SP kiV (Q) - 5P k , v (0)} 



h%k 2 
8tt 



1 



1 



y/SLj/2 + q 2 /4 



CTA 



kd k SP k , v (0) = kd k m% = 



8tt V-cta + fc 2 



f fc 2 
O(k). 



a/-CT A + k< 



0(l/k), 



(18) 



Thus only the mass term m 2 = SP ktV (0) is UV sensitive, while the Q ^ components are not. In the broad resonance 
limit hm — > oo, we therefore find universality from the flow equations: assuming initial conditions 0(1), the loop 
contributions are Oijv 2 ^) and will therefore dominate the physical values of the couplings in the infrared limit, while 
memory of the initial conditions is lost. There is only a single relevant coupling, the mass term. This is the reason 
why we do not have to specify more details of the microscopic inverse boson propagator in Eq. @. For a more 
detailed discussion of universality in the frame of RG equations, we refer to [25]. 

As indicated in Sect. UH the initial condition for the inverse boson propagator is given by Pin iV (0) = v which defines 
the "classical" action for the boson degrees of freedom. Integrating the mass term in Eq. (fTS)) from the initial scale 
fci n — > oo down to the infrared limit k = 0, we find the following relation 



V (B) 



6v in -P , v (0) = ^ k in --^ yf- 



^/(Stt) 



cr A . 



(19) 



The UV renormalization is thus performed by the choice disi n 
Eq. (fT?)) for positive scattering lengths, Pq. v (0) = 0, and conclude the relation (cf. Eq. ([5])) 



Furthermore, we use the exact constraint 



a 



h 2 



1 



(20) 



Since — cta is the gap in the fermion propagator, it may be interpreted as half the binding energy of a molecule, 
£m = 2<ta : — cta is the amount of energy that has to be given to each of the fermions bound in a dimer in order 
to reach the scattering threshold. Therefore, we find the well known relation between scattering length and binding 
energy, 



e M = 2<7 A = -1/Ma 2 . 
The scale dependent inverse boson propagator in vacuum is thus given by 

h 2 



IfcioCQ) = PkAQ) = "(B) + SP Kv (Q) = ^ ( - + \/W2 + g 2 /4-a A + fc 2 ) 



(21) 



(22) 



In the presence of a nonzero binding energy (cta < 0), we may expand the square root. Using Eq. (|20p , we end up 
with 



h 2 v 



(23) 



7 

d k *- — = dk V-*- 



FIG. 1: Graphical representation of the flow equation for the inverse bosonic propagator (cf. Eq. I|17jl) in the vacuum limit. 
Bosons (dimers) are denoted by dashed lines, fermions (atoms) by solid lines. Dots stand for the Yukawa coupling h v , which 
is not renormalized in vacuum. 



which is the inverse Euclidean propagator for elementary bosons of nonrelativistic mass 2M, dressed with a wave 
function renormalization 

h 2 „ h 2 a 

Z v = P= = (24) 

The wave function renormalization coincides with the one obtained in [l7l. [25l[2r3| in the frame of a derivative expansion 
for the effective action. Switching to Minkowski space we obtain the dimer dispersion d> = q 2 / (AM). However, in the 
following calculations we use Eq. (|22|) . 

Let us summarize the results for the two-body problem. We have shown from diagrammatic arguments that the 
fermion propagator and the Feshbach coupling are not renormalized in vacuum. The renormalization of the boson 
propagator can be considered keeping the full frequency and momentum dependence. The flow is driven by fcrmionic 
diagrams only. We emphasize that the solution of the two-body problem is exact as expected for point- like interactions. 
In our formalism, this is reflected by the fact that the two-body sector decouples from the flow equations for the higher 
order vertices: no higher order couplings enter Eq. (117|) . Extending the truncation to even higher order vertices, or 
including the fermion-dimer vertex ipi/jif*^ does not change the situation, since there are only two external lines 
in the two-body problem and the flow equation involves only one-loop diagrams, such that contributions from such 
vertices cannot appear. 



V. THREE-BODY SECTOR: ATOM-DIMER SCATTERING 



In this section we compute the fermion-dimer scattering amplitude A3 from which the fermion-dimer scattering 
length 03 can be extracted in the zero frequency and momentum limit. This problem has been formulated via a 
momentum space inte gral equation a long time ago by Skorniakov and Ter-Martirosian (STM) for more recent 
treatments see @, [13113 Ej. The STM integral equation is derived from a consideration of possible scattering 
processes. These processes form a ladder structure and can thus be resummed in a Lippmann-Schwinger-type self- 
consistency equation, depicted in Fig. [2j 

Here, we present an alternative approach based on a first principles equation, the exact evolution equation for the 
effective action. We derive a flow equation for fermion-dimer scattering, and show how it relates to the STM result. 
Under certain assumptions, we can show the equivalence of both equations. The validity of these assumptions is 
checked via explicit numerical solution of the flow equation. 

The fermion-dimer vertex A3 is computed from 

«.«,,%, <mfh +Q2-Q3- Q t ) = ^^^^ n . ,25, 

In the following we work in the center-of-mass (cm) frame. As our vacuum construction in Sect. IIIII implies, we 
choose the boson to define the zero-point of energy, such that its four-momentum at rest in the cm frame is given 
by P = (0,0). The (Minkowski) cm four-momentum of the fermion at rest reads P = (— a& = — eyi/2, 0), where 
cm is the dimensionless binding energy of the dimer - the fermion propagator is gapped on the BEC side of the 
resonance. Our choice of the zero is different from the one of [2(| where the fermion energy is defined to be zero, 
such that the boson has negative energy cm = 2a a- Of course such a shift in the zero of energy is arbitrary and 
our final equations are independent of this choice. With the cm four-momenta fixed, the effective dependence of the 
fermion-dimer amplitude is reduced to two independent four-momenta, and we will express this fact via the notation 
A 3 (P l5 P 2 ;P). 

We consider the flow equation for the dimensionless fermion-dimer scattering vertex A3 for a specific set of external 
(Minkowski) four-momenta: Consider an incoming fermion with P + P\, an incoming boson with —Pi. The outgoing 
momenta for the scattered fermion and boson can be written as P + P2, — P%. This configuration is in general off-shell 
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3 . As in the STM integral equation, the full off-shell amplitude is needed, since the fermion-dimer vertex also appears 
as a coupling in virtual processes described by one-loop expressions. 

We derive a flow equation for the frequency and momentum dependent fermion-dimer vertex as A3 (Pi, P%\ P). It is 
instructive to consider this equation in a form where the cutoff derivative on the rhs is not yet performed (arising from 
the expansion of the last expression in Eq. ([I])), since this allows for a direct comparison to standard diagrammatic 
techniques, see Fig. El It reads 4 



d k 5X 3 (P u P 2 ;P) = / d k 



1 



5X 3 (P 1 ,Q;P)SX 3 (Q,P 2 ;P) 



[Pf(Q) + Rk lF ) {Pk, v {-Q + p ) + R k lV 

■8X 3 (Q,P2;P)-5X 3 (Pi,Q;P) 



(26) 



K_ 

PA-Q-K 

4 



R 



fe,F 



Pf(~Q ~ P2) + Rk,F 



h 2 



/< 

Q 



p f (-Q -Pi) + Rk,F Pf(-Q - Pi) 
1 



R, 



k,F 



(p F (Q) + R kiF ) (PkA-Q + p ) + Rk >v ) (~M-Q - p i) + Rk,F 



K 



-a 3 (Pi,Q;P) (27) 



P F (-Q - P2) + Rk, F 



We observe the emergence of the running coupling 

X 3 (K U K 2 -P) = 



h 2 



Ppi-Kr - K 2 ) + R ktF 



6X 3 (Q,P 2 ;P) 



-5X 3 {K U K 2 ;P) 



(28) 



on the rhs of the above equation. The rhs thus exhibits a simple quadratic structure in the full vertex A3 = /i^/Pf — SX 3 . 
However, the equation is not a closed equation for A3, since on the lhs only the induced coupling S A3 appears. 

It is possible to make the comparison of Eq. (|26)) to the STM equation directly. However, for the sake of clarity and 
sim ple notation, we first simplify the above equation as appropriate for our purposes. This includes two main steps 
[20I |2l|: (i) Integrating out the loop frequency, and (ii) performing an s-wave projection, since we are only interested 
in low energy scattering here. 

Step (i) is accomplished by noting that only the first fermion propagator Pp(Q) -1 appears with a positive sign for 
the loop momentum Q. For analytically continued (Euclidean) loop frequencies, there is consequently a single pole 
in the upper half plane, while the other propagators (as well as the vertex <5A3, whose pole structure is generated by 
the tree contribution h^/Pp) are analytical in the upper half plane. We may thus perform the frequency integration 
by closing the contour in the upper half plane, implemented by setting the loop frequency u — > — i(q 2 — cta + Rk.F) in 
the remaining momentum space integral. It follows that only values 5A3(Pi = {p\,pi), P 2 = {p^,p 2 ); P) = SX 3 (pi,p 2 ) 
are needed for the solution of Eq. (f2"6"]) . 

At this stage, Eq. (|2"6")l depends on six spatial momentum variables pi,p 2 - However, here we are only interested in 
low energy scattering. For microscopic short-range interactions described by the fermionic scattering length, scattering 
is dominated by the s-wave component [8j. We may account for this fact by performing a suitable projection onto 
the s-wave. For this purpose, we average over the (independent) angles between both incoming and loop momentum, 
and outgoing and loop momentum by applying the integration over 9: 1/2 J , d(cos9). For the inner propagators, 
the averaging is performed explicitly. For the vertex, we define (pi = \pi\) 



SX 3 (p!,p 2 ) 



d(cos9)SX 3 {pi 1 p 2 ). 



(29) 



3 The on-shcll condition reads (u> pi = w P2 , |pi| = |p*2|)- 

4 The flow equation is formulated in Euclidean space, while the physical frequencies are Minkowski frequencies. Therefore, we have to 
analytically continue these Minkowski frequencies ojm to Euclidean frequencies u> and insert these into the flow equation, where we use 
the relation ui = iu>m- 
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With these simplifications, Eq. (f2"T)) reduces to an equation which only depends on the two variables pi,p%. It reads 



d k SX 3 (pi,p 2 ) 



dqq 2 x 
2ir z 



hl(-a-^+R k , lp ^ + 



b 2 



CTA 



Rk,F/2) 



(30) 



x(-5\ 3 (pi,q) + -^-log 



x - 5X 3 (q,p 2 ) + 



log 



Pig + v\ + q 2 - Q"A + Rk ,f 

-pig + p 2 + g 2 - (T A + i?fc,F 
P2<7 + P| + <7 2 - CT A + Rk,F 



dqq 
2 



^{a 3 (pi,9) 



-P29+P2 

8tt 



CTA + -Rfc,F 



hU-a- 



8tt 



4pig 



log 



(7 A + Rk,F 



W - CTA + fc 2 + J4,f/2) 



A3(g,P2) 



-pig + p\ + q 2 - (T A + Rk,F 
8ir\ 3 (pi,q) 



8irX 3 (q,p 2 ) 



p 8tt 
^k.ipj^r 



W - CTA + fc 2 + i?fc, F /2) 



8tt 



|<? 2 - CT A + fe 2 + iZfc, F /2) 



log 



P2<7 + P 2 + <T 



CTA + i?fc,F 



-p 2 q + p 2 



(TA + -R/c.F 



where we have used Eq. (|28|) in the second step, and — <ta = a -2 (cf. Eq. (1201) ). Further we plug in the scale 
dependent solution for the inverse boson propagator obtained from integrating Eq. (j 1 T|) . In the following it will be 
useful to represent Eq. (|30|) in matrix notation 

(31) 
(32) 



d k SX 3 = X 3 -(d k M)-X 3 + (d k L)-M-X 3 + X 3 -M-d k L, X 3 = L — 6X 3 , 
where the matrix elements are denoted by 

h(q, q'), 

L(q,q') 



h 2 

— Mog 

Aqq' -qq 



qq' + q 2 + q' 2 - a A + R kF 



M(q,q') = 



q" - a a 
%Tidqq 2 



R K 



h 2 J-a-i+R k . 



l a 2 

4 y 



(T A + k 2 + R k , F /2) 



The matrix M is diagonal and includes the momentum space measure dqq 2 . Furthermore, matrix multiplication 
means integration over a momentum variable 5 . 

Eq. pop represents a matrix equation which may be solved numerically. However, in the limit where the relative 
cutoff scale c — > oo (cf. Eq. (fTTjO . the equation can be integrated analytically, with a solution that shows the 
equivalence to the STM equation in this limit. 

To see this most straightforwardly, we introduce a rescaled cutoff momentum k' = c^^k and then draw the limit 
c — > oo, kl — const. 6 The bosonic regulator then remains constant, while the implicit scale dependence of the boson 
propagator and the fermionic regulator are suppressed ~ c^ 1 . In particular, for c — > oo, we have 



BvMic- 1 ) -» d k ,M(0), Lie- 1 ) - L(0), Be L^- 1 ) 



0. 



(33) 



For the numerical solution we have discretized Eq. 11311 on a logarithmically spaced two-dimensional momentum grid. 

In order to find /i^-independent results appropriate for the broad resonance limit we need to scale c <x h^/(8n), similar to the detuning 

v. 
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The scale derivative ft reduces to a total derivative in the limit c — > oo. Eq. (|30|) thus reduces to ft<<5A3 = 
A3 • (ft<M) ■ A3, which can be brought in a closed form for A3 since L is fc'-independent in the above limit, 

ft'A 3 = -A 3 • (ft,M) ■ A 3 , or A^ 1 • (ft,A 3 ) • A3 1 = -ft, A3 1 = -ft,M. (34) 

This equation can be integrated straightforwardly with the result 

A 3 = (1+L- M)- 1 -L. (35) 

Here A3 denotes the full fcrmion-dimer vertex in the limit k' — > where all fluctuations are included. Furthermore, 
the flow is initialized with \ 3i k>=k' = £j since for fc' = fc[ n all virtual processes are suppressed and only the tree-level 
graph L is present by construction of the flow equation. 

At this point we observe that quite remarkably, the limit c — > 00 leads to a solution of the flow equation, which is 
independent of the choice of the cutoff function. All regularization scheme dependence drops out in this limit. This 
is because (i) the variables of the differential equation can be separated, and (ii) the function on the rhs of Eq. (|34|) 
becomes a total derivative. This property is lost if c is taken finite. 

We are now in the position to directly compare this result to the STM integral equation [l|, [2(| [2l[ . Expressed in 
the above matrix notation it reads 

\ 3 = L-L-M-\ 3 . (36) 

(For a graphical representation, cf. Fig. \%\) This equation is indeed solved by Eq. (j3"5|) . 

Thus, we reproduce the STM equation in the limit where the relative cutoff scale c — > 00. From the RG point 
of view, this limit implies that the fermions have already been integrated out at the UV scale where the flow is 
initialized. Both the bosonic propagator P Vt k' = Pip,c- 1 k ^ P(p,o an d the tree level graph = L c -i k — > Lq take 
their infrared values in the limit c — > 00 already at the initial stage of the RG evolution. Both these quantities are 
obtained from the elimination of the fermionic degrees of freedom. Hence we recover the STM picture where both 
the boson propagator and the fermionic tree level graph are considered as "fundamental" propagators and couplings, 
in the sense that they serve as the ingredients which allow to construct the Feynman diagrams of an effective theory 
describing fermion-dimer scattering. 

From the RG perspective it appears more natural that fermionic and bosonic degrees of freedom are integrated out 
simultaneously instead of first eliminating the fermions and then the bosons. This is especially true since the bosons 
are not fundamental but build up dynamically as fluctuation induced fermionic bound states. However, the fermions 
are gapped with half the binding energy — <ta = — em/2, while the bosons become massless in the IR limit. In view 
of a judicious choice of c, we may use a simple optimization argument |29j ] . The flow is initialized at k 2 3> —&a- 
In this regime the cutoff function strongly suppresses the flow and the choice of c is obviously arbitrary to a large 
extent. When k 2 > —<ja we observe the onset of nonzero flow. Fermions and bosons couple to each other. The 
optimization sketched below Eq. (jlip requires equal effective mass terms for fermions and dimers k 2 — a a ~ ck 2 . 
Thus c 1 — o~A/k 2 and we conclude a large optimal value of c in the relevant flow regime k 2 < —a a- We emphasize 
that the large value of c is due to the gap —a a in the fermion propagator. This indicates that only small corrections 
can be expected once a full optimization program is implemented. This issue will be addressed in future work. Our 
simple prescription yields 03/0 = 1.22. In Fig. [4] we illustrate the impact of the choice of a fc-independent relative 
cutoff scale on the result for the scattering length ratio a 3 /a for a wide range of c. A c-independent plateau is already 
reached for %-ncjh 2 ^ w 10. The shape of the curve in Fig. 0]in the regime of small c will depend on the choice of the 
cutoff functions, and only in the limit c — > 00 do we obtain regularization scheme independent results. For finite c, 
further optimization - by changing the shape of the cutoff function [2?], [H, [2t| - will reduce the cutoff dependence, 
and presumably also lower the absolute value of 0,3/ a. 

Finally, we specify the relation between the vertex function \ 3 (pi,P2) in the infrared limit k — > 0, and the fermion- 
dimer scattering length a 3 ~ A3(0,0). The proportionality factor is obtained in the following way: First, the boson 
propagator has a nontrivial wave function renormalization Z v = h 2 ,a/(32ir) (cf. Eq. (|24p). Absorbing it into the 
boson field to achieve a standard normalization of the frequency term, we find a renormalized vertex X 3t R = X 3 /Z v . 
The dimensionful version of the vertex (A3 = A3/(2Mfc)) at zero momenta is related to the dimensionful fermion-dimer 
scattering length S3 by the standard relation, \ 3 /Z v — (27r/Af re d) a 3 , where M re d = 2M/3 is the reduced mass of 
atom and dimcr. Expressing these relations in dimcnsionlcss form, we get 



= A 3 , fl (0,0) 

67T 



(37) 
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FIG. 3: Graphical representation of Eq. (|30[) . The shaded circles represent 8X3, the shaded squares A3. The number of the 
corresponding equation in the text is displayed above the equality signs. 




1.0 - 



0.5 - 



0.001 0.01 0.1 1 10 100 1000 

Site/ h 2 v 

FIG. 4: Dependence of the scattering lengths ratio 03/a on the relative cutoff scale c. For c — > 00, the fermions are integrated 
out completely when the flow is initialized and the STM equation (lower dashed line, a-i/a — 1.18) is matched. Furthermore, 
the tree level scattering length ratio 03 /a — 8/3 is indicated (upper dashed line). 

VI. CONCLUSIONS 

In this paper we have studied fermion-dimer scattering in the framework of FRG equations. In a systematic vertex 
expansion for the effective action, we derive the fully momentum dependent flow equations governing the two- and 
three-particle sector of the theory. The sectors form a hierarchy in the sense that the three-particle sector is determined 
by the two-body equations, but does not feed back into the latter. Our flow equations involve diagrams with inner 
fermion lines only as well as mixed diagrams, where both fermion and boson propagators appear in the loop diagrams. 

We investigate in detail the relationship between the FRG equation for the fermion-dimer scattering vertex, and the 
STM integral equation. In our flow equation approach, the starting point is a Yukawa-type theory for fermions and 
auxiliary, non-dynamical bosons after Hubbard-Stratonovic transformation on the level of the classical or microscopic 
theory. We then derive the flow equations for the two- and three particle sector of a theory for fundamental interacting 
fermions from first principles, i.e. a systematic vertex expansion of the exact flow equation for the effective action. In 
the two-particle sector, the dimer propagator is generated gradually in the RG flow. In the three particle sector, we 
encounter a flow equation which is diagrammatically equivalent to the STM equation. In our regularization scheme, 
we can choose the relative scale c at which fermions and bosons are integrated out. Typically, fermionic and bosonic 
degrees of freedom are integrated out simultaneously in RG treatments, accounting for the fact that the dimers are 
not the elementary particles in the original purely fermionic theory. In this work, however, we employ the limit 
c — > 00, where the fermions are integrated out completely already at the initial stage of the RG evolution. Then 
we can integrate the momentum dependent flow equation for the fermion-dimer scattering vertex analytically and 
show the equivalence to the STM equation, implying a ratio of fermion-dimer to fermion-fermion scattering length 
a^/a — 1.18. We argue that the limit of large values for c is physically sensible, since the fermions are gapped by 
half the binding energy, while the dimers become massless for positive scattering lengths (cf. Sect. IIIip . We further 
show that the limit c — > 00 is independent of the choice of the cutoff function. Thus the FRG treatment is suited to 
rederive the STM result correctly. We see this methodolocigal aspect as the main contribution of the present work. 

Furthermore, the techniques developed in this paper can readily be generalized. One direction of future work is 
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to investigate the effects of an effective range, obtained by considering smaller values for the Feshbach coupling h v . 
Furthermore, a fermionic background coupling, corresponding to scattering in the open channel in the atomic physics 
context, may easily be included in our framework. A further goal is the treatment of the four-body problem in 
the presence of a positive scattering length, i.e., dimer-dimer scattering [2(J HH H(| ■ The 1PI four-body scattering 
amplitude {iptp*Lpip*) is described by a dimer-dimer vertex and can be derived along the lines presented above. A 
preliminary analysis reveals that the fermion-dimer vertex feeds into the flow for dimer-dimer scattering. Taking this 
coupling into account produces the same diagrammatic topologies as found in the constructions [20l . l2l| , which go 
beyond the resummation of the bosonic particle-particle ladder performed in (25l . [3l| . Furthermore, the hierarchy of 
the sectors with different particle number is respected: The dimer-dimer coupling (4 particle sector) is affected by the 
couplings in the two- and three particle sectors, but does not feed back in the vacuum limit. 

Our results are also relevant for a quantitatively accurate description of the BEC regime in the BCS-BEC crossover 
problem [24], [H, [HJ : As long as the binding energy of a dimer is the largest energy scale in the problem, the molecules 
should act as point-like entities. Therefore, the many-body physics is expected to be completely determined by the 
value of the dimer-dimer scattering length. As argued above, the dimer-dimer scattering length is in turn affected by 
fermion-dimer scattering. 

Indeed, it has has been demonstrated recently that the value of the dimer-dimer scattering length in vacuum is 
directly reflected in many-body observables at low temperatures [l7j . The formalism provided here is not bound to 
the physical vacuum, and effects of finite temperature and particle density can be straightforwardly included. Our 
technical developments are therefore useful for future studies of the BCS-BEC crossover, and may be seen as the first 
steps towards an approach which combines quantitative precision with a high degree of analytical insight into the 
crossover problem. 
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